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Abstract 

The TASEP {totally asymmetric simple exclusion process) is a basic model for an 
one-dimensional interacting particle system with non-reversible dynamics. Despite 
the simplicity of the model it shows a very rich and interesting behaviour. In this 
paper we study some aspects of the TASEP in discrete time and compare the results 
to the recently obtained results for the TASEP in continuous time. In particular 
we focus on stationary distributions for multi-type models, speeds of second-class 
particles, collision probabilities and the "speed process". In discrete time, jump 
attempts may occur at different sites simultaneously, and the order in which these 
attempts are processed is important; we consider various natural update rules. 
Keywords: TASEP, multi-type, second class particle, speed process. 
AMS 2000 Mathematics Subject Classification: 82C22, 60K35 

1 Introduction 

The TASEP in continuous time was introduced by Spitzer in 1970 (|20)) and can 
be described as follows. It is a Markov process {'r]t)t>o ^^^^ state space E = 
{0, 1}^ where for x G Z we have that site x is occupied with a particle at time t 
iff r]t{x) = 1. Otherwise we say that site x is empty at time t. Starting from some 
initial configuration r]Q £ E, updates occur at each site as a Poisson process of rate 
1, independently; when an update occurs at site x, if there is a particle at site x and 
a hole to its right at site x -|- 1, the particle jumps from site x to site a; -|- 1. If site x 
is empty, or if site x -|- 1 is already occupied, the update has no effect. 

In the model in discrete time, updates occur with some probability /3 £ (0, 1) 
at each site at each time-step. Since updates occur simultaneously, we now have to 
choose an order in which to update the sites. We will consider sequential updates 
(from right to left or from left to right) and sublattice parallel updates (even sites 
first then odd sites). 
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For the model in continuous time there exists a vast amount of hterature. For 
an introduction and background to the topic see Liggett's books |13| (pp. 361-417) 
and |14) (pp. 209-316). However, in some physical models of interest it might be 
more natural to use a discrete time scale. For example in traffic models we can 
consider the reaction time of individuals as a smallest time scale (Blythe and Evans 
[2], Chowdhury, Santen and Schadschneider |3] and Helbing |11| ) and this suggests 
modelling traffic with a model in discrete time. The ASEP (asymmetric simple 
exclusion process, particles jump to the right at rate p and to the left at rate q < p) 
in discrete time was studied for example in Schiitz jl9J, Hinrichsen |12| . Rajewsky, 
Santen, Schadschneider and Schreckenberg |17| and Blythe and Evans [2]. However, 
the behaviour of the models in discrete time has not been analysed in as much depth 
as the model in continuous time. The papers mentioned above are mainly concerned 
with the model on a finite interval with open boundary conditions and just one type 
of particles and analyse density profiles and stationary distributions. 

In this paper we derive further results for the TASEP in discrete time that 
correspond to recently obtained results for the continuous-time model. These include 
stationary distributions for multi-type systems (e.g. [3 [8]), laws of large numbers 
for the path of a second class particle and their connection to competition interfaces 
in competition growth models (e.g. [9l |^), and the TASEP speed process recently 
studied by Amir, Angel and Valko [T]. 

We find that the multi-type invariant distributions for the models with sequential 
updates are identical with those for the model in continuous time, and do not depend 
on the parameter /3. This has the surprising consequence that various collision prob- 
abilities for different particles in a multi-type processes started out of equilibrium, 
of the sort considered in [6] and [T], are also independent of (3 and coincide with 
the values for a continuous-time process. These probabilities correspond to survival 
probabilities of clusters in the associated multi-type competition growth models. At 
the moment, the only argument we have for this property is indirect, using the fact 
that the set of invariant measures is identical for all /?; we do not know of a more 
direct argument based on local dynamics or couplings. 

By contrast, in the case of sublattice-parallel updates, the value of (3 plays an 
important role in the set of stationary distributions. We extend the queue-based 
construction of the multi-type stationary distributions from [3 E] by incorporating 
queues whose arrival and service rates are different at even and odd times. 

The paper is organized as follows. In Section 2 we will give a more formal def- 
inition of the model and introduce the multi-type TASEP. The main results are 
described in Section 3, including results concerning invariant measures and hydro- 
dynamic limits for single-type models which are required in order to state and un- 
derstand the multi-type models described above. The proofs or proof sketches for 
the novel results are found in Section 4. In Section 5 we make some brief remarks 
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about a related discrete-time TASEP model with "fully parallel updates". 



2 Model 

2.1 Models in continuous and discrete time 

The TASEP in continuous time can be described by its generator L. For cylinder 
functions / : — )• M we have 

with the configuration ry^'^+i defined by 

r/(a; + 1) y = x 
r](x) y = X + 1 

Following ideas of Harris (1978) jlOj we can use the following graphical construction 
for the TASEP. Let {(-Pi^)(>o : x G Z} be a family of independent mean 1 Poisson 
processes on a common probability space {Q,A,¥). For x G Z the process marks 
possible jumps from site x: If — Pf_ = 1 and rjt-{x) = 1 then the particle at x 
tries to jump one step to the right at time t. The jump is successful if the adjacent 
site x + 1 was unoccupied, i.e. r]t-{x + l) = 0. Note that for every t > with positive 
probability (e~*) there was no jump in the Poisson process P^ up to time t. Since 
all the Poisson processes are independent there will be infinitely many sites x such 
that there were no jumps in P^. These sites separate Z into intervals of finite length. 
Since no particle can have crossed the boundaries of these intervals, it is enough to 
be able to construct the process separately on each of these finite intervals. 

We can use the same graphical construction to define the TASEP in discrete 
time. All we have to do is replace the family of Poisson processes with a family 
{(-^n)n>o '■ X G Z} of independent Bernoulli processes with parameter (3 £ (0, 1) and 
decide on an update rule for the sites. As mentioned in the introduction we will 
mainly consider the following three update rules: 

• Rule Rl: Updates are processed in order from right to left. 

• Rule R2: Updates are processed in order from left to right. 

• Rule R3: All updates at even sites are processed before all updates at odd 
sites. 

To highlight the difference between the three update rules we can look at the fol- 
lowing example: 
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Figure 1: Configuration at time n and jump marks 
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Figure 2: Configuration at time ri, + 1 if we apply Rl 

Q 
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Figure 3: Configuration at time tt, + 1 if we apply R2 
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Figure 4: Configuration at time n + 1 if we apply R3 



Say we are at time n in the configuration displayed in Figure [H witti particles 
at sites —1 and and holes at sites 1 and 2. There are jump attempts at the sites 
marked with a *. The resulting configurations under the three different update rules 
are as shown in Figures [2] - HI 

Note that in R2, a single particle may jump several times at the same time-step 
(but jumps are only possible onto sites that were already empty at the beginning 
of the time-step). In Rl, several neighbouring particles may jump together at the 
same time-step. There is a natural symmetry between systems Rl and R2 - one is 
transformed into the other by exchanging left and right and exchanging the roles of 
particle and hole. For the last example (R3) the parity of the sites is important. 

In connection with the speed process we will also mention the model with odd /even 
updates (R4). Again this can be obtained from R3 by a simple transformation. 

As seen above, each of these models shows a slightly different behaviour, but if 
we rescale time by a factor f3~^ and let /3 — )• then they converge to the model in 
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continuous time. In this sense the model in discrete time is more general than the 
model in continuous time (which in the following we will denote by RO) since we can 
recover the model in continuous time from the model in discrete time. In discrete 
time we can also consider the model with (fully) parallel updates where all sites are 
updated simultaneously. However, many of the methods developed for the model in 
continuous time that work in the models R1-R3 fail in this case. We will mention 
some questions connected to this model in Section 5. 

2.2 Percolation representations 

Both in continuous and in discrete time, one important feature of the TASEP is its 
connection to last-passage percolation and the corner growth model. Here we consider 
a special case which corresponds to a particular initial condition of the TASEP, in 
which, at time 0, all non-positive sites x < contain a particle and all positive sites 
X > are empty. We label the particles from right to left, so that for i > 1, particle 
i starts at site + 1 at time (and always remains to the right of particle i + 1). 

For n,k > 1, let T(n, k) be the time that particle k jumps to its right for the nth 
time. Then it is well-known that the variables T(n, k) satisfy the recursions 

T{n,k) =mayi{T{n-l,k),T{n,k -1)} + v{n,k) n,k>l (1) 

with boundary conditions T{0,k) = T{n,0) = for all n,k, where v{n,k) are i.i.d. 
exponential random variables with mean 1 . The interpretation is that before particle 
k can make its nth jump, both particle k must have made its (n — l)st jump, and 
particle k—1 must have made its nth jump. Once these two events have happened, an 
amount of time which is exponentially distributed with rate 1 passes before particle 
k makes its nth jump; this is the random variable v{n, k). 

The random variables T(n, k) have an interpretation in terms of last-passage 
percolation times. For an increasing path ir from z € to z' € Z^, i.e. a path 
with increments in {(0, 1), (1,0)}, define the weight of vr by 

S{7r)=Y,v{z"). 

Z"G7T 

Write Il{z, z') for the set of all increasing paths from z to z'; then 

R{z,z')= max 5(7r) (2) 

7ren(z,z') 

is the weight of the heaviest path from z to z' . Then, via the recursions ([T]), it is easy 
to see that T(n, k) = -R((l, 1), (n, k)). In this setting we may interpret the random 
variable f (n, k) as a weight at the lattice point (n, k). 

We turn to the discrete-time case. Now let w{n, k) be i.i.d. random variables 
whose distribution is geometric with parameter /3 G (0, 1) (by which we mean that 
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F[w{z) = k] = {1 - for A; = 0,1,2,...). We define passage-times f{n,k) 
analogous to r(n, k) above by the recursions 

f (n, k) = max |f (n - 1, k), f{n, - 1)| + w{n, k) n,k> 1. 

We will describe three variants on these recursions, which pertain to the different 
update rules Rl, R2 and R3. As above, w{n, k) will correspond to the delay be- 
fore particle k makes its nth jump, once it is free to do so. For i = 1,2,3, let 
T(^(n,k) be the nth jump of particle k under update rule Ki with boundary condi- 
tions r«(0, k) = T^\n, 0) = -1 for all n, k. 

Rule Rl (updates from right to left) 

• Recursions: 

r(^) (n, k) = max jr^^' (n - 1, A;) + 1, T^^) (n. A; - 1) } + w{n, k) 

= f{n, k)+n-l (3) 

• In accordance with the updates from right to left, particles k and A;— 1 can make 
their nth jumps at the same time-step, but two jumps by the same particle 
must be separated by at least one time-step. 

• This corresponds to a percolation model in which as well as weights w{n,k) 
at the vertices (n, /c) G Z^, we have weights of size 1 on each horizontal edge 
between (n — 1, A:) and {n,k). 

Rule R2 (updates from left to right) 

• Recursions: 

r(2)(n, A;) = max {r(2)(n-l,A;),r(2)( n,k - 1) + 1^ +w(n, k) 

= f{n, k) + k-l (4) 

• With updates from left to right, a particle may make several jumps at the same 
time-step, but at least one time-step must separate the nth jump of particles 
A; — 1 and k. 

• In the corresponding percolation model, the weights of size 1 are now on the 
vertical edges of the lattice. 

Rule R3 (even updates then odd updates) 
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Recursions: 



r(3)(n,A;) 



max {r(3)(n -l,k) + l,r(3)(n,fc - 1) + 1} + w{n,k) 

n + k even 

max {r(3)(n- l,/e),T(3)(n,A:- 1)} +u;(n,A;) 

n + k odd 

if{n,k) + ^^±^ n + keven 
= < (5) 
[T{n,k) + ^^^=^ n + A;odd 

• Now the edge weights of size 1 are added to ah edges with an upper/right point 
(n, k) such that n + A; is even. 

For the model in continuous time we have, for x > 0, 

T(\xn],n) , ^ .2 
hm — — ii-^ = (^+l)^ a.s. (6) 

n— >-oo n 

This was essentially first shown in [18] . Replacing the exponential weights by geo- 
metric weights gives 



T{[xn],n) (1 - + 2^(1 - + (1 - /3) 
lim = — a.s.; (7) 

n-5>oo n p 

see for example |16| . Using (I3|)-(l5]), this can easily be used to give similar laws of 
large numbers for T^^\[xn],n), i = 1,2,3. 

We may also view the system as a growth model. For the continuous-time case, 

Gt = {{x,y)e'Ll:T{x,y)<t] (8) 

be the set of vertices whose passage-time is less than t. This gives a cluster in 
which grows over time; there is a 1-1 correspondence relating Gt to the configuration 
of the TASEP at time t\ the length of the row at height k € Z+ is the number of 
jumps particle k has made in the TASEP. In a similar way we can define G^'^\t), 
G^'^\t) and G^^\t) by replacing T in ([8]) by T^^), T^^) or T^^) respectively. 



2.3 Multi-type models 

In the multi-type TASEP each particle belongs to a class y G Z (or more generally 
y € M). All particles can still jump into unoccupied sites. When a particle of class 
k tries to jump into a site which is occupied by a particle of class j two things can 
happen: If A; > j the jump is suppressed and if A; < j then the two particles swap. 
This means that the lower the class of a particle the higher is its priority. 

An A^-type TASEP (containing A^ classes of particles and holes) can be regarded 
as a coupling of A^ ordered single-type TASEPs. If ^q, • • •, f]^ are A^ TASEP 
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configurations such that 7]^ (x) < . . . < tjq (x) for all x € we can use the same 
Poisson or Bernoulli processes (this is called basic coupling) to get a joint realization 
of the TASEPs r]\ ...,r]^. 

The basic coupling preserves the ordering between the processes (since the up- 
dates are processed one by one, this is true for the discrete-time models just as in 
the continuous-time case). Thus we can define a multi-type process ^ by 

N 

6(x)=iV+l- J^7?f(x). 

k=l 

We write = Rvt- Particles of class k occur at sites x where ^(x) = k. For A; > 1, 
these sites represent discrepancies between the processes r]^~^ and t]^. We may 
regard particles of type A'^ + 1 as holes. Then behaves like a multi-type TASEP 
with A^ classes of particles and holes. See for example [7] for further details. 

3 Results 

We will divide this section into three subsections: The first deals with invariant 
measures for single- and multi-type models, the second with hydrodynamic limits 
and the third with multi-type models out of equilibrium. 

3.1 Invariant measures 

Proposition 3.1. For the TASEP in continuous time as well as the discrete time 
TASEPs Rl and R2, the Bernoulli product measures Up with marginals p G [0, 1] are 
the only translation invariant stationary ergodic measures with constant marginals. 
For the TASEP R3, the Bernoulli product measures Hp with marginals p €z [0, 1] on 
even sites and marginals ^i^_p^^ on odd sites are the only stationary ergodic measures 
with marginals that are translation invariant under even shifts. 

Remark 3.1. Interestingly, the marginals of the invariant Bernoulli product measures 
for the models Rl and R2 do not depend on the model parameter /?, and coincide 
with the invariant measures for the model in continuous time. In the model R3 
however, the densities at even and odd sites differ (with a specific relation between 
them) and the measure depends on the parameter /3. 

Proof. For references see for example Liggett [13] for RO, Blytlie and Evans [2] 
for R1,R2 and Rajewsky, Santen, Schadschneider and Schreckenberg [17J for R3. 
The uniqueness statements can be proved following the approach of Mountford and 
Prabhakar [15]. □ 
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We now turn to the construction of invariant measures for systems with more 
than one class of particles. We use the construction based on a system of queues in 
tandem developed in [7j, and begin by recalling notation from that paper. 

Given two processes ai and 02, taking values in {0, 1}^ and representing the ar- 
rival and service processes of a queue respectively, let Z)(ai, Q2) be the process of de- 
partures from the queue. Now define D^^\a) = a, D^'^\ai, 02) = D{ai, 02), and re- 
cursively L>(") (ai, ...,an) = D (L>("-i) (ai, . . . , a-ri—i) ) On) for n > 2. (The process 
can be seen as the departure process from a system of n — 1 queues in tandem). 
Now for a = (ai, . . . ,an) we can define a system of n ordered single-type TASEP 
configurations, denoted Ta = rj = (r/^, . . . , t/") by rj'^ = Z)('^~^+i) (a^, . . . , a„). Then 
the corresponding multi-type configuration = ^(i'---'") is given by ^ = Rr] = RTa, 
with ^(x) = n + 1 — ri^{x) (as in the last paragraph of Section 2). See Remark 

13.21 below for further explanation of the construction. 

We can now state the main result. 

To state this result, we work with systems with jumps from right to left. To return 
to the systems defined before, one simply takes the space-reversal {fjtix) = r]t{—x)). 
Note that time in the queueing system corresponds to space in the particle system. 

Theorem 3.2. If a = (ai, . . . , a^i) ho-s distribution v = Vp^ x ... x Vp^ (fi = 

lip^ X ... X jjLp^ respectively for model R3) with pi < . . . < pn, then the law ofTa = r] 
is invariant for the coupled multi-line TASEPs R0,R1 and R2 (R3 respectively) and 
the law of RTa = Rr] = ^ is invariant for the multi-type TASEPs RO, Rl and 
R2 (R3 respectively) with jumps from right to left. These are the unique stationary 
translation invariant (invariant under even shifts respectively) ergodic measures with 
density p\ of first class particles (density pi of first class particles on even sites), 
density p2 — pi of second class particles ( density p2 — p\ of second class particles on 
even sites), etc. 

Remark 3.2. The mechanism to construct an invariant distribution as described 
above can be depicted in the following way: Take qi as the arrival process and 
a2 as the service process of a queue. Using ai and 02 we can construct a process 
consisting of the departures from this queue (first class particles), unused services 
(second class particles) and times when no service was offered (holes). We then use 
this process as the arrival process for a queue with service process 03 where first 
class particles have priority over second class particles: If there is a service and a 
first and a second class particle are waiting in the queue then the first class particles 
gets served first. In this way we get a resulting process consisting of departures of 
first class particles (first class particles) , departures of second class particles (second 
class particles), unused services (third class particles) and holes. Now we can feed 
this process into a queue with service process 04 and so on. If a = (ai, . . . , a„) has 
distribution v = Vp-^x . . .x Vp^ [p = pp^ x . . . x pp^ respectively) then the distribution 
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Figure 5: Queues in tandem and multi-type configurations ^^^'^^ and ^(^'^'^^ 



of the resuiting multi-type configuration is invariant for the multi-type TASEP. See 
Figure [5] for an illustration. Note that for models RO, Rl, R2, the queues involved 
are simply M/M/1 queues in discrete-time; the same is almost true for R3, except 
that we have different arrival and service rates at odd and even times. 

Remark 3.3. We observe again that the invariant measures for the multi-type TASEPs 
Rl and R2 are the same as the invariant measures for the multi-type TASEP in con- 
tinuous time and that they do not depend on /3. Since the invariant measures for 
the single-type TASEP R3 depend on /3 the same is true for the invariant measures 
for the multi-type TASEP R3. 

3.2 Hydrodynamic limits 

We now move to considering systems out of equilibrium. We consider the particular 
initial configuration given by 

Jl x<0 
[O X > 1 

This "step" initial condition corresponds to the corner growth model and to the 
particular initial conditions for the percolation models described in Section 12.21 We 
define the following functions /q, /i, /2 and /s, which will describe the evolving 
density profile for the continuous-time TASEP and for the TASEPs R1-R3 in discrete 
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time: 



Mu) = l(_oo -ilW + ^(1 --") • 



1-U -T^:/3 



We let as be defined by Mu) = ^ ( a^iu) + M!^KL:_|)^ , so 



asin) = 1 (.) + i (^1 - (2 + n) ^ • M 

For i = 0,1,2,3 let Ti{k,t) (t G IR+ or t G N respectively) be the distribution of 
ivtik + /),/ G Z) in the corresponding model. We have the following result for the 
TASEP in continuous time and the discrete TASEPs R1-R3. 

Theorem 3.3. For any n G R and i = 0,1,2 the measure Ti{[ut],t) converges 
weakly to the Bernoulli product measure with marginals fi{u) and T3{[ut],t) con- 
verges weakly to the Bernoulli product measure with marginals 03(7/) on even sites 
and "^i^J^^u)^^ sites. In particular we have that for any u £ M the limit 

lim(_j>oo IE [?7j(fc)] exists and is equal to fi{u), i = 0,1,2 depending on which model 
we are considering, whenever j tends to u, and liuit^oo^ [vtC^i^])] = o-siu) and 
limt^oo E [r?t(2[|] + 1)] = "l^^a^luffi' ^"^^^ Furthermore, fori = 0, 1, 2, 3, 

the quantities j "^uKkKvt Vt{k) converge a.s. to the constant value fi{w)dw, for 
u < V. 

The first part of the theorem states convergence to local equilibrium: suitably 
rescaled the models converge locally to the unique invariant measures from Theo- 
rem 13.11 This implies the other statements of the theorem. However, in the models 
RO, Rl and R2 we can prove the second part without proving convergence to local 
equilibrium first, while in the model R3 our proof for the second part requires con- 
vergence to local equilibrium. The statements for the model in continuous time were 
proved for the first time by Rost |18) . O'Connell |16| used the connection between 
the TASEP and last-passage percolation to prove an equivalent result about the 
asymptotic shape of the corner growth model (as defined in Section 2). The parts 
of Theorem 13.31 concerning the models in discrete time can be proved using exactly 
the same methods as Rost [18j and O'Connell |16| . In Section 4 we will outline the 
proof for the model R3. 
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Figure 6: Black: Simulation of the corner growth model (R2) for /3 = 0.5 and 
time n = 1000; Red: Limiting shape given by a rescaled version of 
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Remark 3.4. From the convergence to the density profiles /o, /i,/2 and we can 
easily deduce a shape theorem for the corner growth model defined in Section 2. 
The asymptotic shape in the models Rl, R2 and R3 (after rescaling by t) are for 
example given by the functions 



for X E [0,/3], 



52(x) = (v^- V(l-/3)x)' 



(9) 



(10) 



for X G [0, j^], see Figure [6] (simulation with /5 = 0.5 up to time n = 1000), and 



93{^) = (V4x(l-/3)- V4/3-/32x-2/32)^ 



(11) 



for X S 



2/3 
2-/3 



Remark 3.5. We may rescale time as well as space in Theorem 13.31 and look at 
the limit f{u,t) = lim7v-s.oo IE [?7Art(A;)] for limAr_i.oo = u. For the continuous-time 
model, this density profile is governed by Burgers' equation, 



di ^ djfji - /)) ^ 

dt du ' 



(12) 
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the solution, with initial condition /(n, 0) = l(-oo,o] is 

t — u 

f{u,t) = l(_oo + ■ ^[-t,t]iu) 

The differential equation also governs the evolution of the density profile for more 
general initial configurations than the "step" initial condition. We can get equations 
analogous to (jl2p for the models in discrete time. For example, for model Rl, 



h{u, t) = 1 . (n) + i ( 1 - J^^^-^ I • 1 (^) 



solves 



i^r' /? I \ t-u [-i^M 



■5/1, d 0fi{i-fi) 

'W'^ ^ 1-/3/1 



_/l(u,0) = l(_oo,0] (u) 



(13) 



Here ^^I'^p^^^^ = X^n^i /^"/"(l ~ /i) is the probability that a particle jumps from a 
given site to its neighbour in a model in equilibrium with marginal density /i. 



3.3 Multi-type models out of equilibrium 

In this section we consider multi-type TASEPs £ similar to Section 2.3. With 
the results from Theorem 13.31 we can calculate the distribution of the asymptotic 
speed of a single second class particle in the TASEP with initial configuration 



eo(^) 



1 X < -1 

2 X = 

3 X > 1 



As the particles of class 3 are weaker than all other particles in the model we can 
think of these particles as holes. So the second-class particle sees only particles to its 
left and only holes to its right. The second-class particle can be seen as a discrepancy 
between two copies of the "step" initial condition considered in the last section, one 
of which is shifted by one step to the right. Hence the path of the second-class 
particle corresponds to the propagation of the discrepancy under the basic coupling. 
The results for the models in discrete time correspond to the result for the model 
in continuous time first obtained in Ferrari and Kipnis |5]. They prove convergence 
in distribution. In order to prove a.s. convergence we can use the connection to 
last-passage percolation and the growth model: As in Ferrari and Pimentel |9] the 
path of the second class particle corresponds to a competition interface in the growth 
model which has a.s. an asymptotic direction. 
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Theorem 3.4. For i = 0, 1,2,3 let X^^\t) denote the position of the second class 
particle at time t in the corresponding model. Then we have 

t t—^oo 

for random variables with distribution functions 1 — /, for i = 0,1,2 and 03 for 
i = 3. 

Remark 3.6. The proofs for convergence in distribution are analogous to those for the 
model in continuous time (see for example [5]) apart from a small complication in the 
model R3 with the particle-particle coupling. We will give an account of this proof 
in Section 4. For the almost sure convergence we will explain the construction of the 
competition interface for the models in discrete time and prove that the second class 
particle has almost sure an asymptotic speed by using results about semi-infinite 
geodesies in the percolation models similar to [9j. An interesting observation will 
be that the distribution of the asymptotic direction of the competition interface 
corresponding to the path of the second class particle is the same in the models Rl 
and R3, see Remark 14.51 However, this does not imply that the distribution of the 
speed of the second class particle is the same in the two models. 

Remark 3.7. In the continuous model the distribution of the asymptotic speed of 
the second class particle turns out to be uniform on [—1, 1]. The distributions in the 
models in discrete time are more complicated. 

Using Theorem 13.41 we can define the following so-called speed process: Consider 
the multi-type TASEP with initial configuration ^o{n) = n. By Theorem 13.41 we 
know that each particle has a.s. an asymptotic speed Un- Particle n has only stronger 
particles to its left and only weaker particles (that can be seen as holes) to its right 
just like the second class particle in the initial configuration of Theorem 13.41 and 
therefore we can apply Theorem 13.41 to the speed of every particle. We call the 
process {f^n}„gz speed process and denote its distribution by /x. This process is 
stationary, and its marginals (for the various models) are given by the distributions 
in Theorem 13.41 Furthermore, we write Yn{m) instead of im{n) and denote the 
position of particle n at time m by Xn{m) in order to be consistent with the notation 
introduced by Amir, Angel and Valko pj. They have studied the speed process for 
the model in continuous time. Note that both Yn and X„ can be seen as permutations 
of the set Z, and are inverse to each other. We define gi = 1 — fi for i = 0,1,2 
and 53 = 1 — 03, 5(4 = 1 — and our first result is the following theorem 

corresponding to Theorem 1.5 in [J| (note that the labels of particles can now be in 
M instead of just Z): 

Theorem 3.5. For i = 0,1,2, fi^'^\ the distribution of the speed process in model Ri, 
is the unique stationary ergodic measure for the TASEP RO whose marginals have 
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distribution function §{. Correspondingly, fori = 3,4, ^u^*) is the unique stationary 
measure for the TASEP Rj{i), which has marginals distributed according to gi on 
even sites and g^f^^-^ on odd sites, where j(3) = 4 and = 3. 

For i = this gives the result from [Ij saying that the distribution of the speed 
process is itself a stationary ergodic measure for the TASEP in continuous time 
(the marginals are uniform on [—1, 1] in this case). The other parts of Theorem 13.51 
follow from nice dualities between the models Rl and R2 and between the models 
R3 and R4, and the fact that RO, Rl and R2 all have the same set of stationary 
distributions, whatever the value of /3 (as given in Theorem 13. 2p . The dualities are 
given by the following result: 

Theorem 3.6. Consider the starting configuration 1^(0) = n. For i = 0,1,2,3,4 
and any fixed m > the process {Xn\m)^n& has the same distribution as the 
process {Y^^^^\m)}neZ where j(0) = 0, = 2, j{2) = I, i(3) = 4 and = 3. 

The following theorems provide some explicit results about the joint distributions 
of the speeds of adjacent particles (and particles and 2 in model R3). The first 
result is Theorem 1.7 in [1]. The remaining theorems and remarks give analogous 
results for the TASEPs Rl, R2 and R3. 

Theorem 3.7 (TASEP RO). The joint distribution of{Uo, Ui), supported on [-1, if, 
is 

s{x, y)dxdy + r{x)lj^nc=y]dx 

with 

s{x,y) = {l_^ ' and r{x) = ^ 




In particular, P [C/q > Ui] = \,^ [Uq = Ui] = l and P [Uq <Ui] = \. 

Theorem 3.8 (TASEP Rl). The joint distribution of (Uq,Ui) has support on 

2 







and is given by 

si{x, y)dxdy + ri{x)li^.j.=y^dx 



with 

si{x,y) 



i^(l-x) 2(l-y) l=g[(x)g[{y) x>y 



i^(l-x)-i(l-,)-i(Ji^-^i^) x<y 



and 



In particular, P [Uq > Ui] = i, P [Uq = Ui] = ^ and P [Uq < C/i] = i. 
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Remark 3.8. By symmetry, the joint distribution of (C/q, Ui) in the model R2 is the 
same as that of {—Ui, —Uq) in the model Rl. 



Theorem 3.9 (TASEP R3). The joint distribution of {Uo,Ui) has support on 

1 2 

and is given by 



2/3 2/3 
'2-/3' 2-/3 



S2 {x, y)dxdy + r2 {x) l{x=y} dx 



with 



S2ix,y) = < 



X > 



9-i{x)g'i{y) {gA{y) - gA{x)) 

.{2-g^{x)l5-gMP){^) x<y 



and 



r2{x) 



1-/3 f 2(2-/3) 



\ 4-0:2 4 



4^ ) 



In particular, 



>Ui] = ^ (1 -/3)log (y^)) 

1 \ 2(1-/3) 



1-/37 



and 



P[t/o<f/i] 



(l-/3)(2-/3) , 2(1-/3) 



+ 



■log(l-/3) 



^2 ^3 

Remark 3.9. Again by symmetry, we have that under rule R3, ([/i, [/2) has the same 
distribution as {—U\,—Uq). 

Theorem 3.10 (TASEP R3). The joint distribution of {Uo,U2) has support on 

2 



2/3 2/3 
'2-/3' 2-/3 



and is given by 

S3{x, y)dxdy + rz{x)\{^r^=y^dx 



with 



S3ix,y) = < 



9z{x)9z{y) [9a{x) - g^iyY 



2{g4{x)-l){gi{y)^-2g4{x)g4(y)+g4{x)-l 
l-/3g4{x) 

2(g4(y)-l)(g4(x)2-2g4(x)g4(y)+g4(j/)-l _ -,\ 
1-/394(2/) V 



+ 



x>y 



x<y 



and 



„ /™^ _ 93(»)(l-g3W)(l-g4M(l-g4W)) 

2 
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In particular, 

P [Uo >U2] = l 

and 

We see that in every model we have that the speeds are independent on the set 
where Uq > Ui {Uq > U2 respectively). This agrees with the result in continuous 
time. The striking result, shown in [T] for the continuous model, that with positive 
probability the two continuous random variables Uq and Ui are equal, holds also in 
the discrete models. 

Interestingly, the probabilities P [C/q > Ui], ¥[Uo = Ui] and F[Uo < Ui] are the 
same for models RO, Rl and R2, and do not depend on the parameter /3. This is 
rather surprising since /3 is not just a scaling parameter (i.e. we cannot produce 
models with different values of /3 by just applying a time change). In fact, much 
more is true. From the first part of Theorem 13.21 we see that, although the marginal 
distribution of each Ui depends on the model and the value of /3, we can obtain the 
distribution for either of Rl and R2 and any value of /3 by applying an appropriate 
monotone function to each entry Ui (see the proof of Theorem l3.8l for further details). 
Hence the relative ordering of the variables Ui is not affected by the model or the 
value of f3. 

To go further, consider particles i and j with i < j. It's clear that if Ui < Uj then 
particle i can never overtake particle j, while if Ui > Uj then particle i must overtake 
particle j. In [T], it's shown that for the continuous-time model, with probability 1, 
if Ui = Uj then particle i overtakes particle j. The same result can be shown for the 
discrete-time models, although the calculations involved in the argument are rather 
more complicated than those used to prove Theorem 1.14 of [Ij, and we omit them 
here. So, for example, the probability that particle i overtakes particle j is the same 
for models RO, Rl and R2. Indeed, more completely one can define an ordering -< 
on Z by z -< j iff particle j is eventually to the right of particle i. Then we have the 
following result: 

Corollary 3.11. The ordering -< has the same distribution for RO, Rl and R2 and 

for any value of f3. 

It would certainly be interesting to have a more direct understanding of this 
property, based for example on couplings or local dynamics, as well as the indirect 
argument based on the equivalence of multi-type equilibrium distributions. 

Overtaking probabilities in the multi-type TASEP can also be interpreted in 
terms of questions of survival or extinction in multi-type growth models. In [6] a 
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coupling is given between the multi-type TASEP and a three-type version of the 
corner growth model, under which a given cluster survives for ever if and only if 
particle never overtakes particle 1. (The extinction of the cluster occurs if two 
interfaces in the growth model meet ~ the paths of these interfaces are related to 
the paths of the two particles). Different overtaking events in the TASEP can be 
represented by varying the initial condition in the competition growth model. Using 
the results above, we find that the survival probabilities in the growth model will 
remain unchanged if we move from the continuous-time model to natural discrete- 
time models which correspond to models Rl or R2 in the TASEP. Again, this is 
certainly not obvious from the local dynamics of the processes. 

Unlike in models Rl and R2, in the model R3 the probabilities P [f/o > Ui], 
P [Uq = Ui] and P [Uq < Ui] do depend on /3 and the behaviour of the model is 
qualitatively different for different values of f3 (Theorem 13. 9|) : For small /3 we have 
P [Ui <Uo]>F [Ui > Uo], but P [Ui <Uo]<F [Ui > Uq] for large /? (the transition 
occurs at /3 = 0.38860064568 . . .). 

Note however that for /3 — )• 1 the probabilities relating Uq and U2 in Theorem 
13. 101 converge to ^, g and ^, i.e. to the probabilities we get in the continuous model 
and Rl and R2 for the speeds of particle and 1. In a sense, for large /3 the particles 
and 2 in the model R3 behave like adjacent particles in the models Rl and R2. 
This can heuristically be seen in the following way: We consider the particles in 
the model R3 (with large /3 close to 1) starting on even sites. In general, particles 
starting on an even site will move two steps to the right in each time-step since /3 
is large and we update even sites first. If a particle does not jump either during the 
even or the odd update (which happens with probability 2/3(1 — /3)) it ends up on 
an odd site and starts moving left until either 

• (A) it hits a weaker particle to the left by which it cannot be overtaken 

• (B) it does not get jumped over either during an even or odd update because 
the adjacent particle to the left did not try to jump 

In both cases the particle itself will return to an even site (with high probability) 
and resume moving to the right. The particle that caused the stop (either because 
it was weaker or because it did not try to jump) will itself start moving to the left 
until (A) or (B) happens. Now consider the model Rl with large /3. Most particles 
will move one step to the right in each time-step, but some particles do not jump 
and therefore get overtaken until again either (A) or (B) happens (where we remove 
the part "either during and even or odd update"). Particles in these two models 
have different speeds, but the probabilities P [Uq > C/i], P [Uq = f/i], P [Uq < Ui] in 
Rl and P [Uq > ^72], P [Uq = U2\ and P [Uq < U2] in R3 are (almost) the same. 
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4 Proofs 



4.1 Invariant Measures 

As the idea of the proof for Theorem 13.21 is the same for the discrete-time models 
as for the continuous-time model RO we will only sketch the proof. When thinking 
about the model R3 bear in mind that we have different densities on even and odd 
sites. 

Proof sketch for Theorem \3.^ We can proceed in the same way as in [7J : Using ar- 
guments as in |15| we can see that for every parameter p G (0, 1) there exists an 
essentially unique function Hp which maps Bernoulli processes a; on Z x Z onto 
stationary and space-ergodic doubly infinite trajectories {r]n)nez of the TASEP gov- 
erned by uj with time-marginals /Up (Proposition 8 in [7J). For each p and uj we can 
construct a set of dual points Ap(a;) which govern the time-reversal of {r]n)n& and 
again form a Bernoulli process. Also the set of dual points before time m is inde- 
pendent of the configuration r]m (Proposition 10 in [7]). We now take pi < . . . < pn 
and let am = (Om' • • • > '^m) multiline TASEP trajectory governed by to. This 

means that u"' = uj, = Ap^^^ (w'^^^) and (a^)m,ez is the TASEP trajectory gov- 
erned by LJ^ with density p^. Then by the independence of the dual points before 
time m from the configuration r]m we get that the multiline process is stationary with 
product measure ly = x . . . x Vp^ (Proposition 11 in [7]). As in the paragraph pre- 
ceding Theorem 13.21 we define rf = (^rj^, . . . , r/") by r/*^ = (a^, . . . , an)- Then 
induction arguments and some case-by -case checking for n — 2 show that {rj^fiali 
is the TASEP trajectory governed by cj with particle density pk (Proposition 12 in 
[7]) and this implies Theorem 13.21 □ 

Remark 4.1. As mentioned in the beginning of this section, for the model R3 we 
have to think of the pk as densities on even sites and we have to replace the Vp^ by 

Remark 4.2. Inherent in the tandem queue construction for the multi-type stationary 
distribution in model R3 is a version of Burke's theorem for the queues with different 
arrival and service rates on even and odd sites. Consider a queue with arrival process 
An, service process Sn and departure process Let be a Bernoulli process 
with rate pi = (71,72) S (0, 1)^ which means that on even sites arrivals happen with 
probability 71 and on odd sites they happen with probability 72. Motivated by the 
invariant distributions for the TASEP R3 with just one type of particles we want 71 
and 72 to satisfy 

71 (1 - /3) .... 
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where j3 G (0, 1) is the rate at which jumps in the TASEP happen. Analogously, we 
let Sn be a Bernoulli process with rate p2 = 82) S (0, 1)^ where 

(and 7i < 5i, 72 < 82)- The main observation in Burke's Theorem (see for example 
[3]) that shows that arrival and departure process have the same distribution is 
that the queue length process is reversible and that departures look in the reversed 
process like arrivals in the original process. Interestingly, it turns out that in the 
queueing model described above there exists a stationary reversible distribution vr for 
the queue length process which is independent of whether we just observed arrivals 
and services at even sites or at odd sites, vr is given by 

and it is reversible because it satisfies the two systems of detailed balance equations 

vr(j)7. (1 - 6i) = 7T{j + !)(!- 7*) S^ i = 0, 1, . . . 

for i = 1,2 (i = 1 corresponds to even sites, i = 2 corresponds to odd sites). 
This follows from the relations (|14p and ([T5|). As in Burke's Theorem it follows 
from the reversibility of the queue length process that the departure process has 
the same distribution as the arrival process, i.e. D„ is a Bernoulli process with rate 
Pi = (71,72)- 

Indeed, the multi-type construction yields extensions of this result which give 
input-output theorems for the priority queues with more than one type of customer. 
For a discussion of the analogous result in the context of constant arrival and service 
rates, see for example Section 6 of [7]. 

4.2 Hydrodynamic limits 

Proof outline for Theorem \3.3\ for the TASEP R3: The following Propositions cor- 
respond to Propositions 2,3 and 5 in |18| and the proofs are essentially the same as 
in m\. 



Proposition 4.1. For allu (zM. the random variables ^S{[un] , n) = ^ Y2k>[un] Vn{k) 
converge a.s. and in to a constant h^{u), as n goes to infinity. The function /13 is 
decreasing, convex; one has h'i{u) = for u > and h^lu) = —u for u < — 

Proposition 4.2. If h^, is differentiahle at u, one has 



lim E 





'k 




'k' 






2 


) + ??n (2 


2 


-)] 



-2h^{u) 



whenever - tends to u. 
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Proposition 4.3. Let iJL^{k,n) be the distribution of {r]n{k + l),l € Z). If h'^{u) 
exists, any weak limit /ig of the measures /U3 (2 [^] , n) for n — )• 00 is of the form 



with some probability a on [0,1]. Tx is the Bernoulli product measure with density 
b{x) on even sites and density ^^^^^-j^^ on odd sites where b{x) is such that the aver- 
age density is given by x = ^ + ) ■ ^^(^^ means that from Proposition 
4.S\ it follows that the measure a satisfies 



xa{dx) = fsiu) = -h^{u) 



We can use the results from O'Connell |T6] about last-passage percolation (see 
([7])) to calculate the function /13: 



i(2-f -V(4-n2)(l-/3)) 



" - 2-13 



2-/3 



2-/3 



U > 



213 
2-13 



Since /13 is differentiable we can identify fs = —h'^ from Proposition 14.31 as 



fsiu) = -h^{u) 



u < 



1-/3 



2/3 
'2-/3 



.1^ <U< ^ 
2-/3 — " — 2-/3 



U > 



2/3 
2-/3 



As mentioned after Theorem l3.3l in the models RO, Rl and R2 this is enough to prove 
the convergence statements at the end of Theorem 13. 3j this is done using the mono- 
tonicity of the distribution of r]n{k) in k. However, due to the different behaviour 
at odd and even sites, this monotonicity does not hold for the model R3, and the 
average density f^ does not pick up the density fluctuations between even and odd 
sites. Hence without knowing that the model converges to local equilibrium (which 
would allow us to calculate the function 03 from /3) we cannot prove the last part of 
Theorem 13.31 The essential step for proving convergence to local equilibrium is the 
following proposition, the proof of which is again the same as in |18| (Proposition 
6). Let p{k;F;n) =F[r]n{k + i) = l,i e F] for a set F C Z. 

Proposition 4.4. For any finite set F and any e > there exists a 6 > 0, uq such 
that 

/ r7/,m \ / r?lT7,l \ 

< e 





'UW 




un 


;F;n^ 


p(2 


I 2 \ 




2 
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for \u — u\ < 5 and n > uq. Also 

■,F;n + l 





un 




2 





un 










2 





< e 



for < / < [5n] , n > uq. 

Using this proposition and Jensen's inequality we can prove that the measure 
a from Proposition 14.31 is the unit mass on fs^u) and since 6(/3(m)) = as{u) this 
imphes convergence to local equilibrium (see |18] . Section 4 for the details). □ 



4.3 Multi-type models out of equilibrium 

Proof of Theorem \3.4\ ' First we want to outline the proof for convergence of — 



in distribution (as before we have t £ IR+ or t G N depending on the model). This 
follows the ideas in [5j. We want to couple two TASEPs with initial configurations 



1 X < 
X > 1 



and 7?o(^) 




in two different ways and calculate the difference E [^"^([rt], t)] — E [S'^([rt], t)] in 
both couplings. S^{[rt],t) and S'^{[rt,t]) are the number of particles to the right 
of [rt] at time t in r]^ and t]"^. Using basic coupling (i.e. using the same Poisson 
processes {{Pf)t>Q : x G Z} or Bernoulli processes {(fij^)„gi^ : x S Z} for r/^ and r/^) 
gives 



E[S\[rt],t)] -E[S^{[rt],t)] =P xW(t) > [rt] 



(16) 



since we can interpret the discrepancy between r/^ and tj'^ as second class particle 
(see Section 2.3). This works for all models RO, Rl, R2 and R3. The second coupling 
we want to use is called particle-particle coupling. We label the particles in r]Q and 
77q from right to left and let particles with the same label jump at the same time. 
Then under this coupling 

E[SH[rtht)] -E[S\[rt],t)] = P [r?i([rt] + 1) = l] (17) 

in the models RO, Rl and R2. By Theorem 13 . 31 the right hand side of (I17p converges 
to fi{r), so together with (I16p we have 

p[xW(t) > [rt]] > f,{r) 

which proves convergence in distribution for i = 0,1,2. However, in the model R3 
we cannot use the particle-particle coupling as before because this is no longer a real 
coupling. If we let particles with the same label jump at the same time in rj^ and 
rf then the dynamics of the rj^ process are different from the rj^ process: In we 
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update odd sites first and then even sites. We denote by Epp the expectation in rj^ 
and r/^ if particles with the same label jump at the same time in r]^ and ry^ where we 
update in such a way that rj^ is still a TASEP with update rule R3. Then we have 

Epp [SH[rt],t)] =E[SH[rt],t)] 

Notice that starting the second process with updating even sites does not change 
anything as there is no particle on an even site with an adjacent empty site in the 
initial configuration. If we remove the last update of even sites at time t this changes 
the value of 5^([rt],t) if there is a jump from site [rt] to site [rt] + 1 during this 
update. But there can only be a jump from [rt] to [rt] + 1 while updating the even 
sites if [rt] is even. Let us therefore consider odd sites 2[y] + 1 first. We get 



Epp 
and we still have 





'rt 








rt 










= E 


(^2 








2" 









(18) 



Epp 



S 



rt 
~2 



E 



pp 



as before. Hence we get 



X^^\t) > 2 



+ 1 



E 



SU2 



rt 
~2 



m 1 2 



rt 

y 

rt 
~2 



+ l,t 
+ 2 



E 



+ 2 



by the convergence to local equilibrium (Theorem 13.3 
have 



X(3)(t) > 2 



1 



> 



Hence 



lim ] 

t— >oo 



X(3)(t) > 2 

X(3)(t) 



> 



But by monotonicity we 
X(3)(t) > 2 



+ 1 



> r 



Remark 4.3. If the second class particle starts on an odd site (with first class particles 
to the left and holes/third class particles to the right) then 



^ f/(3) 



t—>-oo 



and U^^^ has distribution function ""f^^ ^ accordingly. 
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Figure 7: Pair representation of the second class particle: The figures on the 
left show the system with the pair; the figures on the right show the 
corresponding multi-type system 

Now we want to prove almost sure convergence of the speed of a second class 
particle. Our methods follow the approach in [9]. The idea is to establish a connec- 
tion between the path of the second class particle and a competition interface in the 
corresponding growth model. The cluster in the growth model can be divided into 
two clusters corresponding to events happening to the right and to the left of the 
second class particle. The interface between these two clusters is called the competi- 
tion interface. Using results about semi-infinite geodesies it can be shown that this 
competition interface has almost surely an asymptotic direction. This can be used 
to deduce that the second class particle has almost surely an asymptotic speed and 
since we know the distribution of this speed we can also calculate the distribution of 
the random angle of the competition interface. In the following we will describe this 
method first for the TASEP in continuous time, as given in |9], and then explain the 
adjustments that have to be made in the TASEPs in discrete time. 
In order to establish the connection between the second class particle and the com- 
petition interface we represent the second class particle as a pair consisting of a hole 
and a particle. This reduces our multi-type model to a model consisting only of 
particles and holes and allows us to use the connection to last-passage percolation 
developed in Section 2.2. We let the pair move as follows: If the particle of the 
pair jumps to the right the pair moves to the right (A) and if a particle jumps from 
the left into the hole of the pair the pair moves to the left (B), see Figure [71 Then 
the pair behaves indeed like a second class particle. If we label the particles from 
right to left and the holes from left to right as in Figure [8] then we can consider 
the process {(pn)n>o giving the labels of the pair after the nth jump involving the 
pair. We have ipQ = (1, 1), as initially the pair consists of hole 1 and particle 1, and 
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Figure 8: Pair representation of the second class particle with particles labelled 
from right to left and holes labelled from left to right 



(fn+i — (fn ^ {(0, 1), (1,0)}. (pn Satisfies the recursion formula 

f^n + (1,0) Tiipn + (1, 0)) < r((^„ + (0, 1)) 
[ifn + (0, 1) T{ipn + (1, 0)) > T{ipn + (0, 1)) 



(19) 



If the first label increases the second class particle moves one step to the right and 
if the second label increases the second class particle moves one step to the left. 

Remark 4.4. Note that by inserting an extra site we changed the parity of the sites 
to the right of and including the particle of the pair: The first hole to the right of 
the second class particle is at an odd site while the first hole to the right of the pair 
is at an even site. The second class particle itself is at an even site while the particle 
in the pair is at an odd site. This will be important for model R3. 

Now we want to define geodesies in the last-passage percolation model to define 
a competition interface in the growth model. For z, z' € the heaviest increasing 
path from z to z' (i.e. the path that achieves the maximum in R(z,z')) is called 
the geodesic from z to z' . Note that in the model in continuous time geodesies 
are unique. (In the models in discrete time we will need a rule to break ties to 
achieve uniqueness of the geodesies) . A semi- infinite geodesic starting at z is a path 
TT = (z, zi, Z2, . . .) in such that for every z' = Zj., z" = zi £ n the geodesic from 
z' to z" is (zfc, Zfc+i, 5 • • • 5 -2/) C TT. For a G [0, 90°] a a-geodesic is a semi-infinte 
geodesic with direction a. Now we colour every block Q{i,j) = {i — x (j — 1, j] 
in (R_|_)^\[0, 1]^ either red if the geodesic from (1,1) to passes through (1,2) 

or blue if it passes through (2, 1). The interface between these clusters is called the 
competition interface and an induction argument together with the recursion p9|) 
shows that it is given by the process see Proposition 3 in [9j. Using results about 
the existence and uniqueness of o-geodesics (Propositions 7,8 and 9 in [9]) it can be 
shown that ipn has almost surely an asymptotic direction and we can conclude that 
the second class particle has almost surely an asymptotic speed (Propositions 4 and 
5 in [9]). Now we want to apply these methods to the discrete time TASEPs R1,R2 
and R3. 
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Rl 



The last-passage percolation model for rule Rl was described in Section 12. 2| and in 
particular just after ([3]). Now to adapt to the initial configuration 



r/o(x) 



X < -1 
x = 
X = 1 
x> 2 



we have to remove the '+1' weight from the edge between (1, 1) and (2, 1) and the 
weight from the vertex (1,1) as in the initial configuration we are considering particle 
1 has already jumped over hole 1. We colour Q(l,2) red, Q(2,l) blue and every 
other block Q{i,j) in (M+)^\ [0, 1]^ either red if 

i?«((i,2),(i,i)) >i?«((2,i),(i,j)) 

and blue if 

i?«((i,2),(i,j)) <i?«((2,i),(i,j)) 

(recall the defintion of i? in ([2]); is the corresponding quantity in model Rl with 
the changes mentioned above). This implies that if Q{i,j + 1) is red and Q{i + 1, j) 
is blue, then Q{i + 1, j + 1) is red iff 



r«(i,j + i)>r«(i + i,j) 



(20) 



and blue iff 

f«(i,j + l)<f«(i + l,j) (21) 

where T^^^ is defined as in ([3]) but now in the model Rl with the modifications 
described above. The line (pn^ separating the two clusters is again called the compe- 
tition interface and due to the way we defined the red and blue cluster we have again 
that the competition interface corresponds to the path of the second class particle. 
We can rewrite (|2Up and (|21|) in terms of ipn^ as 



V^^'^ + (i,o) T«((^(,^) + (i,o)) <r«((/7(,^) + (o,i)) 

(^i^) + (0,1) f«(v9« + (1,0)) > f «((/,« + (0,1)) 



(22) 



Note the similarity of (jl9p with ([22]); the difference comes from the fact that in 
the models in discrete time ties are possible. As in [9] we want to prove that this 
competition interface has almost surely an asymptotic direction. First we note that 
the results in [9] about geodesies still hold with geometric weights attached to the 
vertices instead of exponential weights. Secondly, the results still hold in a model 
where ' + 1' weights are attached to every horizontal edge in as these weights do 
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not affect the geodesies (they just give a constant weight to every path from z to 
z', z^z' G The only difference in our case is that there is no weight attached 

to the edge from (1,1) to (2, 1). But this local change does not affect the almost 
sure statements in Propositions 7,8 and 9 in [9]. We conclude that the competition 
interface in our model has almost sure an asymptotic direction and it follows from 
arguments analogous to the ones in [9] that the speed of the second class particle 
converges almost surely. 

R2 

The result for model R2 follows again from the symmetry between Rl and R2. 
R3 

For the purpose of this section it is convenient to change the last-passage percolation 
model corresponding to model R3 a little bit. Instead of updating the even sites 
and then the odd sites during a single time-step, we separate the two batches of 
updates by a half time-step. Then the percolation model corresponding to the initial 
configuration 

Jl x<0 
[0 x>l 

has no weights attached to the edges, while the weights at the vertices are geometric 
with an extra ^ added. As noticed in the beginning of this section, introducing an 
extra site into this model changes the parity of some sites. In order to deal with this 
we attach a single weight to the edge from (1,1) to (1,2) in the percolation 

model corresponding to the initial configuration 

X < -1 
X = 
X = 1 
x>2 

(where we also remove the weight from the vertex (1,1)). This ensures that we apply 

even/odd updates to the left of the particle of the pair and odd/even updates to 

(3) 

the right. Then the movement of the pair ipn corresponds to the movement of the 
second class particle in a model with even/odd updates. 

Again we colour (5(1,2) red, Q{2,1) blue and now every other block Q{i,j) in 
(M+)2\ [0, 1]^ either red if 

^(3)((i,2),(i,j)) >i?(3)((2,l),(i,j)) (23) 
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r/o(x) 



and blue if 

m{{i,2),{i,j))< m{{2,i), (24) 

The interface between these clusters is again the competition interface and is given 
by ifn ■ In terms of (fn we get from (|23p and (j24p that 



V^f) + (1,0) T(3)((^(f) ^ ^^^Q^^ ^ r(3)(c^(f) + (0,1)) 

y.^f) + (0, 1) f (3)(<^(f) ^ ^^^Q^^ ^ f (3)((^(f) + (0, 1)) 



and r(3) are defined as in ([2]) and (l5|) but now we are considering the model 
R3 with the modifications described above. Due to the additional ' + ^' weight on 
the edge from (1,1) to (1,2) the competition interface never encounters any ties in 
this model, i.e. f(3)((^(f) + (1,0)) = f^^^ip'n^ + (0, 1)) does not occur. As in Rl, the 
local change given by the extra ' + ^' edge- weight in this model does not change the 
almost sure statements in Propositions 7,8 and 9 in [9]. The rest of the argument is 
the same as for model Rl. 

□ 

Remark 4.5. We can use the known distributions of the speeds of the second class 
particles (see Theorem 13. 4p together with the hydrodynamic limit results (see The- 
orem [3]3] and Remark 13. 4p to prove the interesting result mentioned in Remark 13.61 
that the distributions of the asymptotic direction of the competition interfaces in 
the models Rl and R3 are the same: 

Proof. The proof exploits the connections made in Proposition 5 in [9]. Similar to 
[9] we let -0?^ = (I^^Ht), J^^Ht)), = {I^'^\t),J^^\t)) be the position of the 

competition interface (i.e. the labels of the pair) at time t and denote by 9^^\ 9^^^ E 
[0, 90°] the random angle of the competition interface for model Rl and R3, i.e. 



4^ 

lim 



(cos(0»),sin(0«) 



for i = 1, 3 



By the arguments in the previous sections we know that these limits exist almost 
surely. Using the asymptotic shape of the growth models given in Remark 13.41 we 
also have 

lim ^ = i»(0«)e^^*'' a.s. for i = 1,3 

t—>-oo t 

where j(*)(0(*)) is the distance from the origin to the intersection of the line given 
by {{u,v) G : tan(0(*)) = ^} and the asymptotic growth interface {x,gi{x)) for 
i = 1,3. With the formulas for gi and in ([9]) and (fTTI) we can calculate j^^\0^^^) 
and j^^\9^^^) explicitly: 



J 



(1) 



/3 



V(l-/3) sin(0(i)) + Vcos(0(i)))' 
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Figure 9: Distribution function for the random angle of the competition inter- 
face in the models Rl and R3 (left) and R2 (right) for (3 = 0.1 (red), 
/? = 0.5 (yellow) and /3 = 0.9 (blue) 



and 



,•(3) 



(0(3)) 



2/3(2-/3) 



(3) Vsm(0(3)) + 2^(1 - (3) cos(e(3))) + /32 cos(0(3)) 



By the connection between the path of the second class particle and the competition 
interface (X»(t) = /«(t) - J^'\t), i = 1,3, since the second class particle moves 
to the right iff the hrst label of the pair increases and to the left iff the second label 
of the pair increases) it follows that 



lim 

t—>-oo 



lim 

t—^oo 



i«(0W) fcos(0(^)) - sin(0«)) = /,(0«) 



t t—^co t 

almost surely for i = 1,3. Using the known distributions of the speed of the second 
class particle in model Rl and R3 (see Theorem 13. 4|) we can calculate 



and 



V^) < a 


= F 


'hiO^'^) > h{a) 


= /i(/i(«)) 


"^(3) < a 


= P 




= 03(^3(0)) 



A calculation shows that 



fiihia)) = asikia)) 



□ 



Proof of Theorem \3.6l Let the operators (T„ be defined by 

(TnY Yn<Yn+l 
(TnY = < 

[ Y otherwise 

where Tn exchanges Yn and Yn+i in Y . The proof of the following general Lemma is 
the same as in [1] (Lemma 3.1). 
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Lemma 4.5. For a fixed sequence ii, . . . , in Z we have 

d , ^_l 

We have the relation Yx„(m)('^) = between the X and the y process. Since 
/? < 1 each site has a positive probabihty that no jump occurs at that site at any 
given time. At each time-step these sites separate Z into finite intervals and the 
events on these intervals during that time-step are independent. In the model with 
updates from right to left we apply a finite sequence of operators Oi^ - • • ai^. where 
ii, . . . ,ik is an increasing sequence (since we update from right to left). Lemma [4.51 
states that applying cjj^ • • • cjjj. is the same (in distribution) as applying cjj^ • • • 
(i.e. updating from left to right) and taking the inverse permutation. But given the 
configuration 1^(0) and performing the updates from left to right we get 1^ (1) 

(2) 

and this is the inverse permutation of Xn (1). So 

Inductively we get that this holds for all m > 0. The other three parts follow in 
the same way. In the case with even/odd updates, ii, . . . ,ik is a sequence such that 
there exists al<j<A; + l such that ii is odd for / < j and ii is even for I > j. □ 

In order to prove Theorem 13.51 we will need the following Lemma. We state it 
here for the model Rl, but analogous results hold for the other models as well. The 
Lemma corresponds to Lemma 4.1 in [Ij. 

Lemma 4.6. Consider two TASEPs, Y^^^ and Y^^\ as functions of the same 
Bernoulli points on Z x N (i.e. under basic coupling). We set y^^^(O) = n and 
1^^"'^'*(0) = aj ■ ■ ■ (Tj^kYn^\o) for some j G Z and k > 0. Let {Un^} be the speed 
process ofY^^^ and {Un^} be the speed process ofY^^K Then = aj^k ■ ■ ■ 

Proof. Every other particle than {j, . . . j -\-k + l} is either stronger than all particles 
{j, . . . , j + k + 1} or weaker than all particles {j, . . . ,j + k + 1}. Any swap of a 
particle other than {j, . . . , j + k + 1} will happen in both Y^^^ and Y^^\ So for 
any i ^ {j, . . . , j + k + 1} we have X^^\m) = X^^\m) for all m > and therefore 
U^^^ = U^^^ for those i. In Y^^'^ particle j + k -\- 1 is always to the left of all other 
particles {j, ...,j + k}. So = min{C/j^^ . . . t^j+fe+i}- Define j < r < j + k + l 

by 

min{c/f)...^«,^J = ^« 
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Then for i = r,r + 1, . . . ,j -\- k we have = U-^_^\ and for i = j, j + 1, . . . , r — 1 



Proof of Theorem \3.5l Consider a BernouUi process on Z x Z. Half of this process 
{7j X N) is used to construct the TASEP Y^^\ For any Z G Z we can translate 
the Bernoulli process by / (i.e. take points of the form (n, m + I) where (n, m) is 
in the original process). We can restrict this translated process to Z x N and use 
this restricted process to construct another TASEP. Let U'^^\l) = {u!;^\l)} be the 
speed process for the TASEP that has been constructed using the Bernoulli process 
translated by /. For every /, U^^\l) has distribution So we have to show that 

{Uri\l)^ behaves like a TASEP with updates from left to right. In order to do 
this we look at a transition {Un\l)^ {Un\l + 1)}. The effect on the original 
TASEP of changing from translating by / to translating by / + 1 is that some finite 
sequences of a operators of the form (Jj ■ ■ ■ Uj+fc are added to be applied to the 
TASEP before the original sequence of operations. At each location a a operator 
is added with probability (3. The previous Theorem shows, that applying each of 
these finite sequences has the same effect on the speeds as applying each sequence in 
reverse order to the speed process. This shows that {u'f' {I)} behaves like a TASEP 
with updates from left to right and therefore the measure /x^^^ is stationary for the 
TASEP with updates from left to right. 

The proofs for the other three models are essentially the same (using the appropriate 
versions of Lemma 14. 6p . □ 

The following Lemma will allow us to do the explicit calculations for the joint 
densities of the speeds in Theorems 13.81 - [3. 101 using the connection between queueing 
models and the invariant measures introduced in Theorem 13. 2i Here D{gi) is the 
domain of the distribution function gi (i = 0, 1, 2, 3, 4). 

Lemma 4.7. If F : D{gi) — )• {!,..., A^} is non- decreasing then for the TASEP 
n^m the distribution of {-F(f/n^)}n is the unique ergodic stationary mea- 
sure of the multi-type TASEP model Rj{i) with types {1, . . . , A^} and densities A/ = 
c/j(sup{F"^(/)}) - 5j(inf{F"i(/)}) for type l = l,...,N (j as m Theorem\3^. 

Proof. The proof is analogous to the proof of Corollary 5.4 in pLj. □ 




This shows that U^-^^ 



O'j+k 



□ 
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With the help of this Lemma we can do ah the calculations needed for the results 
in Theorems 13.81 - 13.101 Depending on the model we are considering we will choose 
the function F from Lemma 14.71 to be Fi = mm{j : gi{u) < Xj} for some increasing 
sequence (xi, . . . , XAr_i) in [0,1]. Then we put Vn = F{Un) and the distribution 
of the Vs is given by the invariant measure for the multi-type models and can be 
calculated explicitly using the queueing representation. 

Proof of Theorem \!J.8l By Lemma r4.7l we have with = 3 that Fi (Un) is distributed 
according to the unique ergodic stationary measure of a 3-type TASEP with updates 
from left to right and densities 

Ai = P [xq < gi{Un) < xi] = xi - Xq = Xi 
A2 = P [xi < gi{Un) < X2] = X2- Xi 

A3 = P [X2 < gi{Un) < X3] = X3 - X2 = 1 - X2 

Ai is the density of first class particles, A2 is the density of second class particles and 
A3 is the density of third class particles (or holes). Recall that Vn = Fi[Un)- Using 
the queueing representation for the unique ergodic stationary measure of a 3-type 
TASEP we can calculate the joint distribution (Vq, Vi) explicitly. This distribution 
depends on the Xj. Taking suitable derivatives with respect to these Xj we get the 
density of the corresponding speeds. We have for example 

P [C/o < g^\xi) <Ui< gi\x2)] =f[Vo = 1, Vi = 2] 

= X1X2 (X2 - Xi) 

since the probability of having a second class particle at position 1 is X2 — xi (since 
this is the density of second class particles) and to have a first class particle at 
position we then have to have an arrival (probability xi) and a service (probability 
X2) because having a second class particle at site 1 means that the queue was empty 
at that time (so in order to have a departure at site we need an arrival at site 0). 
Remember that in Theorem 13.21 the particles in the TASEP jumped from the right 
to the left. If we want to consider the TASEP with jumps from the left to the right 
(and that is what we are doing here) we have to read the queues from right to left. 
So position 1 comes before position and the probability of having a second class 
particle at position 1 is independent of arrivals and services at position 0. So for 
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uq < Ul we put xi = gi{uo) and X2 = 9i{ui) and get as density 

dxi dx2 d d 
P (7o G duo,Ui e dui\ = - — - — - — - — xiX2 {x2 - xi) 

duQ du\ dx\ dx2 

= (1 - ^^o)"^ (1 - (251 («i) - 251 (uo)) 

1-/3 3 3 

= ^(i-norMi-^in 



1 — txi V 1 ~ '"0 



Similarly, we have 



P [5r'(xi) < C/i < 9l^{x2) < i7o] = P [Vb = 3, = 2] 

= (1 - X2) (X2 - Xi) 

and therefore we get as density for uq > u\ (putting xi = g\{u\), X2 = giiuo)) 
¥[Uoeduo,U,edu,]=(-p-) (-p.)-^-^{i-x2){x2-x,) 



duo J \ dui J dxi dx2 
= ^(l-^o)-^(l-ui)-^ 
= 9'i(.uo)g[{ui) 

To get the density for uo = ui we consider 

P [g^\xi) < Uo, Ul < gTHx2)] = P [Vb = 2, Vi = 2] 

= (1 - Xi)x2 (X2 - Xi) 

and let xi,X2 — )■ giiu). We get 

F[Uo,Uiedu]= lim 



xi,X2^gi{u) ^{X2) — ffi ^(xi) 

vr^(i-gi(n))gi(n) 

2/3(l-u)t 



2^2(i_^)lV 2/32(l-t.fV/3 
2/33 (l-n)^ 

To get the probabilities in the Theorem we only have to integrate the densities over 
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the appropriate ranges of uq, ui and u. Alternatively we can use the following: 



(^)l =p[5i(^,S'^)<5i(C/J'^)" 

= P [[/(°) < C/f 

_ 1 

~ 3 

(*) follows from the fact that the distribution of {gi{UrP)} is the unique translation 
invariant stationary ergodic measure for the TASEP R2 with marginals uniform on 
[0, 1]. The distribution of {go{Un^)} is the unique translation invariant stationary er- 
godic measure for the TASEP in continuous time. Since the stationary distributions 
for the multi-type TASEPs RO and R2 are the same (see Theorem [321) {OiiUn^)} 
has the same distribution as {go{Un^)}. □ 



The proofs for Theorems 13.91 - 13.101 work in exactly the same way. 



5 Fully parallel updates 

Finally we mention the model with "fully parallel updates". If an update occurs at 
site X at time t (which happens with probability f3 as usual), this update causes 
a jump from x to x + 1 if and only if r]t-i{x) = 1 and rjt-i{x + 1) = (that is, 
the jump is already possible before any other updates at the current time-step are 
performed) . 

There are several important differences between this model and the models we 
have studied earlier. The Bernoulli product measures Vp are no longer invariant. 
Furthermore, the basic coupling no longer preserves an ordering between different 
initial configurations, and so it is no longer clear how to define a multi-class system. 
If we use basic coupling to couple two systems which start with one discrepancy 
at the origin then this single discrepancy can generate additional discrepancies. It 
would already be interesting to know how the leftmost and rightmost discrepancies 
behave. Do they have asymptotic speeds, and if so are the speeds random or deter- 
ministic? There is still a natural percolation representation, and we can still obtain 
a hydrodynamic limit result in the sense that ^Ylun<k<vnVt{''^) converges a.s. to 
the constant value f{w)dw, for u < v and some function /, but the stronger 
result that lim„_>oo E [rjnik]] exists and is equal to f{u) whenever ^ tends to u does 
not follow using the same methods as in the proof of Theorem 13.31 
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